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I. INTRODUCTION 



A few years ago, a RCNP group including two of the present authors , has 
started to systematically study the properties of instanton systems with respect to 
confinement by using multi-instanton simulations. The instanton is a quasi-pointlike 
topological object in Euclidean space-time and, being selfdual, satisfies the Yang-Mills 
field equations || . Mixed instanton-antiinstanton configurations || (gas or liquid) only 
locally are approximate solutions, but they offer the attractive possibility to solve, in 
a semiclassical-like fashion, the £7a(1) problem (explain the large rf mass) || and 
to explain chiral symmetry breaking (See Ref. || for a recent review.) One 

main focus of our previous work on instantons, neglected by others, was to explore 
under which circumstances a dense, however random instanton system could provide 
confinement. In the beginning ]I] it has been pointed out that the monopoles induced in 
the instanton system and detected by Abelian projection have the clustering property 
held necessary for confinement. As for the QQ force itself, in Ref. [0,|3| quantitatively 
the conclusion has been reached, after some refinements, that with the widely accepted 
density of 1 fm~ of instantons plus antiinstantons and an average instanton size of 
p = 0.4 fm, only 40 % of the static QQ force could be reproduced at distances around 
R « 1 fm. 

This result was achieved using an instanton size distribution based on the idea 
of freezing of the strong coupling constant at large distance. This assumption leads to 
a dimensionally dictated behavior of the size distribution like dn(p)/dp oc p~ 5 in the 
infrared. One should remember that the main parameter of any such model, the in- 
stanton density and the instanton size distribution are phenomenological input, beyond 
justification coming from a truly semiclassical approximation. The latter is available 
only for a single instanton ||10|| , however afflicted with the famous infrared divergence. 

Phenomenologically, in the instanton liquid model instantons are occupying Eu- 
clidean space with a density of n = N/V = 1 fm -4 , and with an average radius 
p — 1/3 fm one gets a packing fraction / = n p 4 »i (1/3) 4 . For our purposes we fix 
only an average instanton radius p = 0.4 fm. We study the influence of the instanton 
density, considering 1 fm -4 as the "physical" value. This leaves some room for the 
choice of the size distribution. 

Lattice search for instantons is the method of choice to obtain these quantities 
and the size distribution relying on first principles. As for pure gluodynamics, only 
lattice studies are at hand presently to quantify the instanton structure of the vacuum. 
Over the last years the results have not much converged (compare the recent conference 
reports Ref. [I I - jllfl ). Different groups roughly agree on the size of instantons within a 



factor of two, e.g. p = 0.3. ..0.6 fm for SU(3). There is no agreement at all concerning 
the density N/V which is strongly dependent on the procedure to remove size 0(a) 
vacuum fluctuations (cooling, smoothing etc.) and technical parameters (cooling steps, 
cooling radius etc.)Q As a tendency, lattice studies give higher density and larger 



* Sometimes the density N/V = 1 fm 4 is used as a criterium to stop cooling, in this way 

2 



instantons than phenomenologically assumed. It is fair to say that only the topological 
i 

susceptibility, Xtop/V& = 0.45(3) for SU(3) and 0.50(2) for SU(2), and the average 
instanton size are supported by lattice studies. The instanton size distribution is much 
more delicate to assess. We will comment on that later. 

This paper is dedicated to a critical reappraisal of the instanton model to de- 
scribe, besides other features of the Yang-Mills vacuum, its confinement property, and 
we shall clearly point out the deficits of this model. Although seemingly insufficient 
in quantitative respect, it might be interesting to see to what extent the contribution 
of instantons to the confining force depends on the p distribution. Also the question 
whether the instanton generated QQ force can be reproduced in accordance to Abelian, 
monopole or vortex dominance needs some clarification. One has to answer the ques- 
tion at which scale the model can be replaced by an effectively Abelian model with 
condensing Abelian defects. Monopole and Z(N) vortex mechanisms are presently the 
leading candidates for an effective infrared description comprising confinement. We 
shall see that these descriptions are applicable to the instanton mediated force as well. 
This is expected, corroborating earlier work jnj and the widely-studied interrelation 



between instantons and monopoles [|T^,|Ty] and the newer studies concerning the inter- 
play of instantons and vortices |17|,|18| , respectively. 



Said in another way, what we want to clarify is the complementarity between 
the explicit semiclassical-like description in terms of continuous instanton fields on 
one hand and the monopole and vortex aspects on the other. The latter degrees of 
freedom seem to become physically dominant in the infrared. This is the place where, 
for our purposes, the lattice discretization appears: it is an infrared matching scale 
between the instanton picture and the gauge singularities which become manifest in 
the result of gauge fixing and Abelian projection. If the discretization scale is chosen 
too small, monopole and vortex degrees of freedom can be identified as well. However, a 
complementary description of the instanton mediated force can be achieved only if the 
matching scale is somewhere between the size of and the distance between instantons. 

For the sake of clarity we stress that this paper is not a lattice study. The lattice 
is employed here only to enable the necessary coarse-graining of a continuum model. 
In a lattice gauge theory investigation the role of the discretization scale a would be 
completely different. There a is an ultraviolet cut-off which permits, at the cost of a 
running bare lattice coupling (3 oc 1/g 2 , not to deal explicitly with fluctuations of smaller 
and smaller wave lengths. In this case, the requirement of scale invariance in the limit 
aA oc exp ^jfP \ — ^ is indeed crucial and must be confirmed for any dimensionful 
quantity to be physical. Instantons can be found which are stochastically generated 
within the sample of fields. Being lumps of gauge invariant topological density they are 
of immediate physical importance. As mentioned above, their average size is relatively 
well-defined, independently of the lattice scale a (for suitable methods of suppressing 
the shortest fluctuations living on the lattice). For the confining Abelian defects the 
situation is somewhat different. They are identified by gauge fixing (which is practically 



specifying the remaining instanton characteristics. 
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performed on the lattice of scale a), and the scale invariance of the corresponding 
density, of the distribution of length or area etc. are controversial. Moreover, it is 
generally agreed that in order make these defects condensing they have to be defined 
with some extension (blocked monopoles or thick vortices), a scale which becomes 
decoupled, in the continuum limit, from the lattice scale a. As far as instantons are 
discussed as a possible microscopic mechanism to induce Abelian defects, finding the 
correct matching scale is tantamount to define this extension. 



There have been also statements in the literature |19[ criticizing the instanton 
contribution for violating Casimir scaling of the QQ force. While this feature of the 
non-perturbative force at intermediate distances |30[ is largely not understood within 



Yang-Mills theory and continues to pose a problem for other models of confinement, 
we want to make clear to what extent this criticism is justified in the case of instanton 
based models. Finally, many discussions on the instanton generated QQ force are based 
on a dilute gas formula worked out by Callan, Dashen and Gross ||21|| . Therefore it 



is of interest to demonstrate under which circumstances the result of multi-instanton 
simulations deviates from the one-instanton approximations made in Ref. ||21[| . In the 
physical range of packings the effect of different instantons on the Wilson loop is not 
expected to factorize anymore. 

We believe that this study contains some lessons for other attempts of an explicit 
modeling of the QCD vacuum. The observed scale sensitivity of the monopole or center 
vortex description seems to be a more general feature of the complementarity between 
semiclassical continuum models for non-perturbative vacuum structure and condensing 
defects. It should be remembered that the concept of Abelian and monopole dominance 
had been introduced as a property of gauge fields in the infrared p2|-pT| , not necessarily 
on the lattice. Practically, however, all evidence comes from doing the gauge fixing and 
Abelian projection on the lattice for gauge fields generated by lattice simulation. In the 
present work the instanton model continuum configurations will be discretized with the 
purpose to perform the same steps following pB],|2B[] for the monopole part and [ESH3T] 



for the vortex part of the heavy quark force. Explicitly calculating these contributions, 
it turns out that the discretization scale a is an influential infrared coarse-graining 
parameter and must be chosen in correspondence to the typical size and density of the 
disordering continuum non-perturbative configurations forming the vacuum. This is 
what our instanton model clearly illustrates. 

In section 2 we describe how the lattice discretization of continuum instanton- 
antiinstanton configurations is done. In section 3 we demonstrate that the one- 
instanton description of the QQ force breaks down already in moderately dense in- 
stanton gases. In the same section we study to what extent the forces for quarks and 
adjoint charges deviate from Casimir scaling. Section 4 deals with the question how 
coarse-graining the lattice has to be chosen such that Abelian projection, the monopole 
gas and the center projection are able to reproduce the non- Abelian instanton-mediated 
force. In section 5 we show, for the case of monopoles, that this necessary blocking 
corresponds to the percolation property held necessary for confinement. We summarize 
and conclude in section 6. 
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II. CONTINUUM MULTI-INSTANTON CONFIGURATIONS 
AND THEIR DISCRETIZATION 



We base our studies of the QQ force and possible complementary descriptions 
in terms of monopoles and center vortices on a model which comprises the Yang-Mills 
vacuum as an ensemble of random collections of instantons and antiinstantons. The 
interaction is partly taken into account in the size distribution given below. 

In order to fix the scale of our continuum model, we choose an average instanton 
size p = 0.4 fm. The average size is most solidly defined by the profile of the topological 
density after a few cooling steps. Our choice realistically applies to SU(2) Yang-Mills 
theory. Unless stated otherwise (when we study the density dependence of the instanton 
mediated force in section 3) the density is chosen N/V = (Nj + Nj)/V = 1 fm~ 4 with 
the packing fraction / = p 4 y = 0.0256. 



We adopt the sum ansatz [p2 



AM = J2 A i( x ~i z k, Pk, o k ) + J2 A l( x 'i z h Ph Ok), (i) 

k k 

in terms of instanton and anti-instanton solutions in the singular gauge where an 
instanton is written as 

Here p and z denote the instanton size and the space-time position of the instanton 
center, respectively. The instanton solution can be rotated in color space by the color 
orientation matrix O. The 't Hooft symbol ff^ is defined as ff^ u = e b ' J ' u (l — 5 AfJ- 5 4u ) — 
6 bn 6 4» + §bv 5 4» The anti . selfdual so i u tion A\ is obtained replacing ff^ by if^ = 

(— iy" 4 + sv4 fj b v. The instanton solutions have several gluonic collective modes related 
to variations of parameters like size and position (five collective degrees of freedom). For 
pure SU(2) gauge theory, the color orientation matrix is characterized by 3 parameters 
(the Euler angles). For N c colors the number of collective parameters (and gluonic zero 
modes) generalizes to 4N C . 

Actually, we generate the ensemble of instantons and anti-instantons by ran- 
domly placing Zk in a 4-dimensional Euclidean continuum box. The (adjoint) color 
orientations Ok are taken randomly with the Haar measure, and the instanton sizes pk 
are sampled according to the following size distribution : 

dn(p) , 



dp 



a/- 5 exp(-/3p7p2) ( 3 ) 



with b = lliVe/3. Here a and (3 are fixed by normalizing to the space-time density as 
Jq°° dn(p) = N/V and the average size / °° p dn(p) = p N/V. In the explicit configura- 
tions, the instanton number Nj is taken equal to that of the antiinstantons Nj. 
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This form of the size distribution was originally established in a mean-field treat- 
ment of color-averaged interactions ||32|| . Averaged over orientations, the interaction of 
equal-sign and opposite-sign pairs is repulsive, S int oc p\ p\. In the mean-field approx- 
imation for the one-instanton distribution, this results in a suppression exponentially 



in p 2 . 



As for the lattice evidence for such an instanton size distribution we refer to the 
systematic problems mentioned in the Introduction. One should keep in mind that 
instanton finding algorithms depend on cooling or smoothing. Huge instantons distin- 
guished by a very weak field strength could become washed out within the first steps 



of cooling together with the noise. Smith and Teper 03] have reported an instan- 
ton p distribution from quenched SU(3) lattice ensembles resembling an exponential 
cutoff at large p. The size distribution found by Hasenfratz and Meter |R3 has been 



interpreted by Shuryak [|35| in an exponentially suppressed parameterization in terms 



of a dual Higgs mechanism. There are also other warnings, motivated by the wealth 
of semiclassical configurations known today ||37|| , not too early to take for granted the 



size cut-off as seen in the lattice studies of Refs. 33 



Very recently, in Ref. P8| , a distribution like (0) has been demonstrated to 
emerge, in dimension d > 2, from a very general model of soft (inflatable) spheres with 
excluded- volume interaction which inherits from QCD not more than the semiclassical 
perturbative instanton fugacity proportional to p b ~ 5 . In the grand canonical Monte 
Carlo simulation of this system the p dependent fugacity has been used in conjunction 
with the excluded- volume constraint. 

Because the instantons are randomly placed inside the box, it might be interest- 
ing to know how frequently instantons are strongly overlapping, which we define by 
pi p2/Af 2 > 1 (where A i2 is the Euclidean distance). For the size distribution (|3|) and 
the physical density of 1 fm -4 , the number N< of strongly overlapping pairs (including 
all possible pairings between instantons and antiinstantons) relative to the total num- 
ber N = N + + N_ of instantons and antiinstantons amounts to N</N = 0.072. This 
fraction actually depends on the width of the size distribution. For instance, for a sharp 
distribution ^j^- oc 5(p — p) one would find, with the same density and p, a smaller 
fraction N < /N = 0.056 of closely packed pairs. We are aware of the problem that the 
linear superposition ansatz should be replaced by some better ansatz (ratio ansatz) for 
close instanton-instanton pairs and by the streamline ansatz for instanton-antiinstanton 
pairs at high densities. For the observed percentage of close pairs in the physical range 
of densities, we have neglected this potential complication. This improvement would 
become more important with an algebraically falling size distribution. 

In our actual calculation, we cover the random multi-instanton configuration by 
a lattice similar to Ref. ||. For the calculation of the fully non-Abelian force the 
discretization actually would not be necessary. In fact, we choose a sufficiently fine 
lattice spacing of a = 0.05 fm (this should be compared with the average instanton size 
of p = 0.4 fm). For this case (of the finest lattice) the link variables are given simply 
as U^(x) = exp[iav4 /J (x + \pa)] using the gluon field of the multi-instanton system on 
the mid-point of the link I = {x, p}. In connection with Abelian gauges and Abelian 
projections, we will apply a coarser discretization scale. Then the SU(2) link variables 
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are constructed by integrating the vector potential defined on the continuum space as 

fx+fia 



U = UJx) 



P exp 



/ 



P n exp [ a A^(x + f 10 ) 

3=1 



(4) 



Here, the path from x to x+pa has been subdivided into smaller segments with step size 
of a = a/ 1, and the above path ordered exponential has been calculated as a product 
over finer links defined on these segments. Actually, we take a = 0.05 fm as the segment 
size. This construction becomes increasingly important when the discretization scale 
of the lattice gets comparable with or larger than the average instanton size. 

Periodicity of the lattice gauge field configurations has been enforced by placing 
the 4-dimensional Euclidean box of size V = (4.8 fm) 4 (to be covered by the lattice) into 
a bigger box. The 3-dimensional boundary in each of the eight directions is extended to 
a 4-dimensional slab of thickness 1.2 fm continuing the basic box. In each of these slabs 
copies (phantom instantons) are placed of instantons which are near to the opposite 
boundary. These are also included into the sum (|l[) representing the continuum vector 
potential. Then, along the links restricted to the basic box the above construction is 
performed. 



III. THE QQ POTENTIAL FROM THE RANDOM INSTANTON LIQUID 

AND THE DILUTE GAS LIMIT 



The effect of instantons on the heavy quark potential has been first discussed 
in Ref. |H] by using the dilute gas approximation. There a formula has been derived 
according to which the static potential could be written as 



V(R) 



dp VP, 



(5) 



where dn(p)/dp is the respective instanton size distribution and the factor two refers 
to instantons plus antiinstantons. Here, a dimensionless potential W(R/p) (R = \r\) 
is defined by integration over the instanton position x 

W (-) = — \ / d 3 x tr [ U(r-x)U\-x) - 1 ] . (6) 
\P J 2 P 6 J 

The fundamental static charges (j = |) are represented by static trajectories positioned 
at and r in 3-space (with R = | r |). The color matrix U(f—x) represents the path 
ordered exponential along one of the infinite trajectories traversing the instanton and is 
given as U (r*— x) = exp[i n r-(f—x) / [ (f—x) 2 +p 2 ] 1 ^ 2 ]. It depends only on the instanton 
size p and the shortest distance from the trajectory to the instanton center. Notice that 
the instanton orientation does not matter. That part of integration is normalized to 
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unity. Replacing the trace |tr(C/i Ul) in (|) by the trace in the adjoint representation, 

|tr a d(C/i U\) = \ ((tr(Ui U\)) 2 — lj one can extend @ to the adjoint (j = 1) charges. 
For the comparison in this section we adopt the instanton size distribution in the form 
of eq. (|J) and confront the dilute gas result (|5]) with multi-instanton simulation results 
where the same distribution has been used for sampling of the p^'s. In this section, the 
total density N/V is a free parameter. 

The potential (for fundamental and adjoint charges) obtained from eq. ([5]) is 
simply proportional to the density N/V as shown in Fig. [T] by the solid lines. The 
data symbols in Fig. p] denote the results of an uncorrelated multi-instanton simula- 
tions. The comparison is made for the instanton densities (a) N/V = 0.05 far , (b) 
N/V = 0.20 fnT 4 , (c) N/V = 0.60 far 4 , and (d) N/V = 1.00 far 4 . The results 
of the simulation are based on a statistics of (a) 100, (b) 300, (c) 500, and (d) 1000 
configurations respectively. The left panels show the potential between fundamental 
charges, while the right panels between adjoint charges. 

The simulation data have been obtained, using the lattice discretization of the 
random instanton-antiinstanton configurations as described above, from expectation 
values of Wilson loops 

(W full (R,T)) = (Tr Yl^) (7) 

zee 



where the contour C is a rectangular closed path of size R x T on the finest lattice of 
a = 0.05 fm, a spacing almost 10 times smaller than the average instanton size. The 
non-Abelian potential 

V(*) = -Ilog. < H/ '»»< fl ' r » 



(w M (ii,r-i)> 



has been constructed from the discrete logarithmic time derivative. 

We see that eq. (|J) is only approximate and applicable only for a system of 
widely separated uncorrelated instantons and/or antiinstantons. Then the effect of 
single instantons exponentiates in a Wilson loop at finite temporal extension T. Only 
at very small distance, where overlapping of instantons can be neglected, simulation 
and analytical formula agree in the quadratic behavior of the potential. 

Already for a moderate density of N/V = 0.2 fm -4 deviations are remarkable 
which develop differently with higher density for fundamental and adjoint charges. 
From the comparison we conclude that the dilute gas formula for fundamental and 
adjoint charges is valid only up to a density of 0.05 fm -4 which corresponds to a 
packing fraction f — 1.3 x 10~ 3 . Corrections become visible first at large distances of 
R w 1.5 fm. We show this for the density 0.2 fm -4 corresponding to a packing fraction 
/ = 5.2 x 10~ 3 . The resulting potential from the multi-instanton simulation continues 
to rise where the dilute-gas formula starts already leveling off. 

For the comparison of the instanton contributions to the forces between " quark" 
and "antiquark" in the fundamental and adjoint representation, respectively, it is more 
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illuminating to consider the derivatives of the potentials shown in Fig. |I[ We are 
interested now in a instanton gas of realistic density and show in Fig. the result of 
simulations for N/V = 0.6 fm -4 only, together with the dilute gas result. The ratio 
between the adjoint and the fundamental charge is < 2 at all distances and is closest 
to the Casimir ratio, | at R — 0.5 fm « p where the instanton contribution itself is 
maximal. The ratio becomes quite small as the distance increases. This tendency is 
seen for the other densities, too. 

In Fig. [| we show the density dependence of the QQ force at two distances, 
R = 0.9 fm « 2 p and R — 1.5 fm, above (a) for fundamental, below (b) for adjoint 
charges. It is interesting to see for the smaller distance that both forces rise with 
the instanton density N/V almost linearly, with a slope compatible with the dilute 
gas (one-instanton) approximation, up to N/V ~ 0.6 fm -4 . The forces drop down 
significantly at higher density, as shown by simulations at N/V = 1 fm -4 , in the case 
of fundamental charges relatively to the dilute gas extrapolation, while in the case of 
adjoint charges even absolutely. A similar behavior is obtained for all R < 1 fm. On 
the other hand, for R > 1 fm the behavior is different. We show this for R ~ 1.5 fm. 
The force in the fundamental representation rises almost linearly in the density with a 
slope exceeding the dilute gas approximation by rs 50 %, while the force in the adjoint 
representation drops already at small density below the dilute gas result. 

IV. VARIOUS ABELIAN PROJECTIONS 
OF MULTI-INSTANTON CONFIGURATIONS 

In this section, we want to clarify how concepts like Abelian dominance, monopole 
dominance, center dominance can be applied to semiclassical-like multi-instanton con- 
figurations. We shall compare the fully non-Abelian force discussed in section 2 with 
the results obtained after standard techniques of gauge fixing and Abelian projection 
have been applied to the multi-instanton fields and the respective contributions to the 
static QQ force have been evaluated. Technically, as a coarse-graining device, the con- 
tinuum configurations are latticized, and the results will depend critically on the lattice 
spacing compared with the average instanton size. 

For this study the periodic Euclidean box of size V = (4.8 fm) 4 is discretized 
with four different lattice spacings, a = 0.2 fm, a = 0.4 fm, a = 0.6 fm, and a = 0.8 fm, 
corresponding to lattices 24 4 , 12 4 , 8 4 and 6 4 , respectively. In all cases, we have chosen 
a = 0.05 fm as the segment size (integration step) to construct the links (path ordered 
exponential) of the lattice configuration. How periodicity is enforced has been described 
in section 2. 

First, we consider the maximally Abelian gauge [[Mj (MAG), which is defined by 
maximizing the functional 

Rma = E tr P^x) r 3 U^x) r 3 ] (9) 
with U^(x) = U°(x) +iT t U t (x). In the MA gauge, the SU(2) link variable U^x) is 
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decomposed as 




M^(x)u^(x) 



1- | c M (x) | 2 -c*Jx) \ ( e iB ^ 



c M (x) J I- I c^z) 





(10) 



where the Abelian angle variable O^x) and the non-Abelian variable c^(x) are defined 
in terms of U^x) as tan^(x) = U%(x)/U%(x), c^(x)e ie ^ = [-U*(x) + iUfa)]. To 
clarify the contribution of Abelian components to the static force, we consider the 
Abelian projection of full non-Abelian link variables to the Abelian ones u^. This is 
tantamount to replacing, in a Yang-Mills vacuum configuration put into MA gauge, of 
the SU(2) link variables by U(l) link variables, U^{x) —> u^{x). Before this is done the 
off-diagonal parts UHx) and U^(x) of gluon fields have been minimized by the gauge 
transformation which has to be found iteratively. 

One can further decompose the diagonal gluon component 9^ into the monopole 
part 9™° and the photon part 9^ h . Using the forward derivative d fl f(x) = f(x + fia) — 
f(x), the 2-form 9 t±l/ (x) = d fl 6 u (x)—d u 9 fl (x) defines the field strength which is separated 
as follows 

e iiV (x) = e ia/ {x)+2Tm ia/ (x) (11) 

into a gauge invariant regular field strength O^ix) = mod2 7r 6 l MV (x) G (— 7r,7r] and a 
Dirac string part n^ u (x) G Z. 

From each part of the field strength the photon part 9^ h (x) and the monopole 
part 9™°(x) of the U(l) vector potential can be reconstructed, for instance 

9™(x) = 2tt£ n-\x - x')d v n^(x'), (12) 



using the lattice Coulomb propagator □ 1 = (d li d') 1 , where d'^ denotes the backward 
derivative. Then, we can construct the Abelian projected Wilson loop as 



{W abel (C)} = (e*p[if O^x) dxp\) (13) 

which contains, besides of a "photonic" Wilson loop the monopole projected Wilson 
loop 

(W mon (C)) = (exp[z j£ 9™ dx,\) (14) 

uniquely defined factor. 

On the other hand, we can also consider the center projection of an instanton 
configuration. This requires to go through the lattice discretization, too. Then, we 
use the direct maximal center gauge [p9|1 . This gauge is defined as the gauge which 
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brings directly the full link variables U\ as close as possible to the center elements ±7 
by maximizing 

Rmc = E ( tr U,{x)} 2 . (15) 

X,fM 

In this gauge the center projection is defined by 

Z^x) = sign[tr U^x)]. (16) 

After the gauge fixing, there remains only a local Z{2) symmetry as U^x) — > z(x)z(x + 
fia)U^(x) with z(x) G Z{2) Then, we can construct the center projected Wilson loop 

as 

(W center (C)) = (l[Z l ). (17) 

l£C 

In Fig. |], we show the instanton mediated non-perturbative force for the stan- 
dard density 1 fm -4 and fixed average instanton size p = 0.4 fm. Instanton sizes 
were sampled according to the exponentially damped size distribution (0). The sim- 
ulation data are based on a statistics of 1000 multiinstanton configurations for each 
discretization scale. 

As shown in Fig. f|(a), for a lattice spacing much smaller than the typical size 
of instantons, the Abelian force calculated after Abelian projection from flUf ) does not 
reproduce the non- Abelian force. Moreover, the Abelian force is practically reproduced 
after the monopole component of the Abelian field has been removed (by the "photonic" 
Wilson loops alone). The monopole contribution to the heavy charge force turns out 
to be smaller by a factor of two. The Z(2) force calculated within the center-projected 
configurations is completely negligible. 

For comparison we show in Fig. [|(b) the same for a coarser lattice of lattice 
spacing a = 0.4 fm = p. In this case, the Abelian force has increased but it is still 
far from reproducing the non-Abelian force. The center component of the force has 
become measurable and is comparable with the monopole component which has also 
increased somewhat. 

This trend continues when we consider a even coarser lattice in Fig. ||(c) with 
a = 0.6 fm = 1.5p. In Fig. @(d) we show the coarsest lattice with a = 0.8 fm = 2p. 
Now the slopes of the potentials calculated in the various projections are in agreement 
with Abelian, monopole and center dominance, with an ordering of the quasi-string 
tensions osu(2) > CAbciian > o"z(2) ~ cTmono- Here, of course, the potential can be looked 
at only at one or two lattice spacings. Therefore the values of the string tension should 
be considered only as rough estimates. 

V. MONOPOLE CURRENT PERCOLATION 

The monopole current can be defined using u^x) following DeGrand and Tou- 
ssaint J40|. Since the Bianchi identity regarding the Abelian field strength 9^ u (x) is 
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broken by the decomposition in eq. (|TT|), a monopole current k^x) can be defined on 
the dual link {*x, p} as 



k^'x) = —e^ a pd v e a p(x + fia) = -d, 



where h^ u (x) = \s ^ ua pnap(x + pa). This current is obviously conserved, d^k^(*s) = 
which results in closed monopole loops on the 4-dimensional dual lattice. Here, & 
denotes the backward derivative. 

In the previous paper [0], "block-spin" transformations (of type II) have been 
applied to the monopole currents [f4T| , f42"|1 in order to study how the global structure of 



the monopole loops changes for Abelian projected multi-instanton configurations. For 
large-scale blocking, leading to monopole currents defined on a very coarse (infrared) 
lattice, monopole clustering appears, i.e. clusters of big monopole length are seen to 
form. 

Instead of this, in the present study, we consider first the Abelian projection 
performed on lattices with different lattice spacings. The different lattices would be 
related, in the sense of by type I block spin transformations. For us, the difference is 
just in the choice of the discretization scale for our multiinstanton configurations. The 
histograms have been, however, obtained in independent runs. In the case of Figs. || (a) 
and (b), 2000 multiinstanton configurations have been evaluated, in the case of Figs. 
[I (c) and (d) the statistics was 500 configurations each. The Figs. [5] show histograms 
dN/dl monopole clusters normalized to the average occurrence in one configuration. 
The integral over lengths, J dN(l), gives the average number of separate clusters per 
configuration, while J I dN(l) gives the average total length. When the projection is 
performed on a lattice with a << p (in our case a = ~ p) the histogram of monopole 
loops per configuration with respect to their loop length (cluster size) contains only a 
component of relatively short loops as shown in Fig. |§(a). If the projection is done on 
a lattice with a = 2p the histogram Fig. |^(c) has already a component of monopole 
currents which percolate. One very long monopole cluster of complicated structure 
appears per configuration with a discretization scale a = 2p, similar to what has been 
observed previously ||, as shown in Fig. |5|(d). 

In Fig. ^| we show, in a 3-dimensional cut, the instanton and monopole content 
of a characteristic configuration at that density. The cross-shaped endcaps of the bars 
symbolize that the monopole world lines leave (enter) the volume going to (coming 
from) another timeslice. In MAG, the monopole network (percolating at the coarsest 
lattice) is wrapping the instantons. 



VI. SUMMARY AND CONCLUSIONS 



We have studied the static quark potential induced by a random instanton liquid. 
We have made the present study in the sequel of the previous work || where an infrared 
suppression of the instanton size distribution dn(p) oc p~ v dp has been assumed. In 
this paper, we adopt a size distribution as dn(p)/dp = ap b ~ 5 exp(— [3p 2 / p 2 ) (suggested 
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by lattice instanton searches) adapted to a space-time density of N/V = 1 fm -4 and 
an average instanton size p = 0.4 fm. We find that the gross effects do not depend on 
the detailed form of the size distribution. 

First, we have compared the non-Abelian potential inferred from the Wilson 
loop embedded in the simulated random multi-instanton ensemble with the well-known 
dilute-gas expression which reduces the potential to the effect of isolated single instan- 
tons. We found that the dilute gas expression differs essentially from the simulation 
results if the density is bigger than i of the phenomenological value, due to many 
instantons interacting with a Wilson loop in a non-exponentiable way. This explains 
why the quasi-linear part of the QQ force at R ~ 1 — 1.5 fm needs to be found from 
a simulation. We have found that the force between adjoint charges behaves more 
complicated in the multi-instanton system. 

Second, we have investigated the Abelian and center projection of the gauge 
field configurations built up by random instanton configuration. The gauge fixing 
and subsequent projection has been performed introducing lattices of different lattice 
spacing, and we have studied the effect of changing the discretization scale, in particular 
checking lattice spacings a = 0.2, 0.4, 0.6 and 0.8 fm. In the finer lattice case, the 
monopole contribution to the heavy charge force turns out to be much smaller than the 
Abelian projected force. The force calculated for the center-projected configurations is 
completely negligible. With larger and larger discretization scale, Abelian, monopole 
and center dominance becomes restored, which is almost perfectly illustrated for a = 2p. 
Then the quasi-linear non-Abelian force can be almost reproduced by the corresponding 
projected Wilson loops. 

We stress this result, although the random instanton liquid model as such does 
not appear to be a realistic model for the confining aspect of the Yang-Mills vacuum, 
because this seems to reflect a more general feature to be kept in mind for more 
sophisticated semiclassical models. 

We see that in the multi-instanton liquid, at a = 2 p, Abelian dominance of 
the heavy charge force amounts to about 90 %, which is mainly due to the Abelian 
monopole component. The singular part of the Abelian gauge potential accounts for 
80 % of the Abelian force. The center projected force now also describes about 80 % 
of the non-Abelian one. 

We have considered the monopole loop percolation behavior in the random multi- 
instanton configurations using the same typical lattice spacings. Only after performing 
the projection with a discretization scale a > p we could find a percolating component 
among the monopole currents, represented by one or few clusters of huge loop length 
per configuration. This parallels the previous observation in Ref. || that blocking of 
monopole currents is necessary to see percolation/condensation happening. 
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VII. FIGURE CAPTION 
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FIG. 1. The static non-Abelian potential derived from simulations of the random in- 
stanton liquid and from a dilute gas (one-instanton) approximation for fundamental (i) and 
adjoint (ii) charges The comparison is made for instanton densities N/V = 0.05, 0.2, 0.6 
and 1.0 fm -4 . The curves describe the dilute gas approximation, the data points simulation 
results. 
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FIG. 2. The static non-Abelian force derived from simulations of the random instanton 
liquid and from a dilute gas (one-instanton) approximation for fundamental (circle) and 
adjoint (square) charges. The comparison is made for an instanton density N/V = 0.6 fm - . 
The curves and the dotted curves describe the dilute gas approximation for fundamental and 
adjoint charges, respectively. 
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FIG. 3. The density dependence of the forces at distance R = 0.9 fm and R = 1.5 fm for 
(a) fundamental and (b) adjoint charges. The solid curves describe the dilute gas approxi- 
mation and the data points correspond to the simulation results. 
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FIG. 4. The non-Abelian static potential of fundamental charges compared with the 
Abelian projection, its monopole component and with the potential in Z{2) center projection, 
for different lattice spacings (a) a = 0.20 fm, (b) a = 0.40 fm, (c) a = 0.60 fm and (d) 
a = 0.80 fm. The instanton size and density is fixed to p = 0.4 fm and N/V = 1 fm. 
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FIG. 5. Histograms of connected monopole clusters per configuration with respect to 
their lengths shown for different discretization scales (a) a = 0.20 fm (b) a = 0.40 fm, (c) 
a = 0.60 fm and (d) a = 0.80 fm. 
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